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Abstract. This paper is a direct continuation of pQ. The Hamiltonian aspects of the theory are 
further developed. Within the framework provided by the first paper, the problem of minimality for 
constrained calculus of variations is analyzed among the class of differentiable curves. A necessary 
and sufficient condition for minimality is proved. 
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Introduction 

This is the earlier of two papers concerned with a further analysis of the arguments of pQ on the 
topic of constrained calculus of variations. 

Consider an abstract system 23, subject to a set of differentiable and possibly non-holonomic 
constraints. In its essential features, the problem tackled in these papers is the one of finding, 
within the class of the admissible evolutions 7 of 23 joining two given configurations, the one which 
minimizes the action functional 

J[ 7 ] := [j?(t,q\t),z A (t))dt 

being Jzf (t, q 1 , . . . , q n j Z ; . . . j Z ) £1 function on the admissible velocity space of the system and 7 
the lift of 7 to the latter. 

A preliminary step in this direction was just provided by [I], in which the stationarity conditions 
for the functional I [7] were studied through the analysis of its first variation. The algorithm resulted 
into a set of conditions characterizing stationary curves among the class of piecewise differentiable 
ones, i.e. among the class of continuous curves having a finite number of discontinuities in their first 
derivative. Furthermore, the analysis pointed out the equivalence between the given constrained 
Lagrangian problem and a free Hamiltonian one. 

The present paper is instead focussed on establishing whether a given extremal curve provides 
a local minimum of the functional I [7] . For the time being, the argument is developed within the 
only class of differentiable curves, leaving to the latter paper all problems arising from the possible 
presence of corners. 

As far as the present work is concerned, the demonstrative strategy may be briefly sketched as 
below. 

On the first instance, in §[21 taking a given extremal curve 7 into account, a sufficient condition 
for minimality is proved by means of the analysis of the second variation of the action functional. 
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In this connection, the geometrical tools developed in pQ prove very useful in order to simplify 
the matter. The argument turns out to be subject to the solvability along the whole of 7 of a 
Riccati-like differential equation. Unfortunately, because of the non-linearity of the latter, Cauchy 
theorem can only guarantee the existence of purely local solutions. Hence, at this stage, the global 
solvability is leaved among the hypotheses. 

Subsequently, in £j3l the attention is shifted onto Jacobi vector fields. They represent a special 
class of infinitesimal deformations which link families of extremal curves. They are used to inves- 
tigate the processes of focalization and, by means of the further concept of conjugate point, to give 
a necessary condition for minimality. 

At last, in Jfl the sufficient condition and the necessary one are merged together, showing that the 
lack of conjugate points along the extremal curve 7 implies the solvability of the above non-linear 
differential equation along the whole of it. 



1.1. Preliminaries. This Section is meant to be a reference tool consisting of a brief review of 
those contents of [1] that will be involved in the subsequent discussion. All results are stated 
without any proofs nor comments. Of course, the reader is referred to pQ and references therein for 
a thorough description of the subject. 

(i) Let V„+i R denote a fibre bundle over the real line, henceforth called the event space, and 
referred to local fibered coordinates t, g 1 , . . . ,q n . 

Every section 7 : R — > V n +i will be interpreted as the evolution of an abstract system OS, 
parameterized in terms of the independent variable t. The first jet-bundle ji(V n +i) — > V n +i, 
referred to local jet-coordinates t,q l ,q l , will be called the velocity space. The first jet-extension of 
7 will be denoted by .71(7) : R — > ji(V n +i) ■ 

We shall indicate by V(V n +i) — > V n +i the vertical bundle over V n +i — > R and by 
V*(y n+ i) V n +i the associated dual bundle, commonly referred to as the phase space. By 
definition, the latter is canonically isomorphic to the quotient of the cotangent bundle T*(V n +i) by 
the equivalence relation 



To make it easier, we preserve the notation ( , } for the pairing between V(V n +i) and V*(V n+ i) . In 

this way, every local coordinate system t, q l in V n +i induces fibered coordinates t, q % ,p% in V*(V n +i) , 
uniquely defined by the requirement 



1. Geometric setup 






(ii) The pull-back of the phase space through the map ji(V n +i) 
commutative diagram 



TT 



described by the 



C(ji(V n+ i)) 



* V*(V n+1 ) 



(1.1) 



c 



7T 



il(V n +l) 



7T 



* V n+ i 
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gives rise to a manifold C(ji(V n +i)), henceforth referred to as the contact bundle. 

We shall refer it to fibered coordinates t, q l , q l ,Pi , related in an obvious way to those in ji(V n +i) 
and in V*(V n +i) • Every a € C(ji(V n +i)) will be called a contact 1-form over ji(V n +i). 
Notice that, by construction, C(ji(V n +i)) is at the same time a vector bundle over ji(V n +l), with 
projection £, and an affine bundle over V*(V n +i), with projection k. 

The bundle C(ji(V n +i)) — > Ji(V n +i) is canonically isomorphic to the vector subbundle of the 
cotangent space T*(ji(V n +i)) locally spanned by the 1-forms 

(1.2) a/ := dq i - tfdt i = l,...,n 

With this identification, for every a € C(ji(V n +i)) the coordinates Pi(o~) coincide with the 
components involved in the representation 

(1.3) a = pj(cr) oj\c,(,a) 

Eventually, it is worth recalling that the contact bundle carries a distinguished Liouville 1-form 
© , locally expressed as 

(1.4) Q = Pi J = Pi {dq l -q l dt) 

(iii) The presence of differentiable constraints restricting the admissible evolutions of the system 
is formally accounted for by a commutative diagram of the form 

A — Jl(Vn+l) 

(1.5a) 

V n +1 ^=^= Vn+1 

where 

• A ——> V n +i is a fiber bundle, representing the totality of admissible velocities; 

• the map A ji(V n +i) is an imbedding; 

• a section 7 : R — > V n +i is admissible if and only if its jet-extension ^1(7) factors through A, 
i.e. if and only if there exists a section 7 : R — > A satisfying j'i(7) = i ■ 7. Under the stated 
circumstance the section 7, commonly referred to as the lift of 7, is called an admissible 
section of A. 

Referring the submanifold A to fibered local coordinates t, q\ . . . , q n , z\ . . . , z r , the imbedding 
i : A — > ii(V n +i) is thus locally represented as 

(1.5b) q i = ^ i (t,q 1 ,...,q n ,z\...,z r ) i = l,...,n 

while the admissibility conditions for a section 7 : q % = q l (t), z A = z A (t) read 

(1.6) ^ = ^{t,q\t),...,q n (t),z\t),...,z r (t)) 



^For simplicity, the same notation is used for the 1-forms H1.2|) as well as for their pull-back on C(ji(V n +i)) ■ 
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Every section a : V n +i — > A is called a coniro/@ for the system. The composite map i ■ a is called 
an admissible velocity field. In local coordinates we have the representations 



(1.7a) 
(1.7b) 



a : z A = z A (t,q 1 , ...,q n ) 
i-a : q* = r(t,q\z A (t,q*)) 



A section 7 : II — > V n +i and a control cr : V n +i — > *4 are said to belong to each other if and only 
if the composition a ■ 7 coincides with the lift 7 : It — »■ A- 

(iv) The concepts of vertical vector and contact 1-form are easily adapted to the submanifold A : 
as usual, the vertical bundle V(A) is the kernel of the push-forward T(A) — ^ T(V n +i) while the 
contact bundle C(A) is the pull-back on A of the contact bundle C(ji(V n +i)), as expressed by the 
commutative diagram 



c(A) — U ccn(v n +i)) 



(i. 



y4 



-> Jl(Vn+l) 



The manifolds V(A), C(A) will be respectively referred to coordinates t,q l ,z A ,v A and t,q l ,z A pi 
In this way, setting 



(1.9) 

we have the representations 



Cj 1 := i*{uj l ) = dq l - ip l (t, q k , z A ) dt 



X = v A (X) 



d 

dz A 



Vie V(A) 



(7 



Pi(a) uj 1 | C(ct) V<TGC(i) 



The pull-back @ := i*(Q) of the 1-form (jl.4j) will be called the Liouville 1-form of C(A). In 
coordinates, eqs. (|1.4p . (|1.9p provide the representation 



(1.10) 



Q = Pico 1 = pi (dq l - ifj 1 dt) 



Finally, setting straightforward composition of diagrams (|1.8|) , (| 1 . 1 j) yields the bundle 

morphism 



(1.11) 



C(^) — ^ F*(V n+ i) 
C 

^ — V„+i 



expressing the contact bundle C(^4) as a fibre bundle over V*(V n +i), identical to the pull-back of 
the latter through the map A — > V n +\ ■ 



The term control is intuitively clear since every assignment of the section a determines the evolution of the 
system from a given initial configuration through the solution of a corresponding Cauchy problem for the ordinary 
differential equations (|1.6|) . 
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1.2. Vector bundles along sections. From here on we shall stick on the geometrical environment 
based on diagram (|1.8p . 

(i) Given any admissible section 7 : El — > V n +i lifting to a section 7 : R — > A, we denote by 
V("f) — > R the vector bundle formed by the totality of vertical vectors along 7 and by -A (7) — > R 
the analogous bundle formed by the totality of vectors along 7 annihilating the 1-form dt. 

Both bundles will be referred to fibered coordinates — respectively t, u % in V(j) and t, u\ v A in 
A (7) — according to the prescriptions 

(1.12a) <=► X = n i (X)(A') 

(1.12b) XeA( 7 ) <=► * = _ + . 

Recalling that the first jet-extension of the section 7 is denoted by .71(7) , we have the following 
Proposition 1.1. The first jet-bundle jiiV^)) is canonically isomorphic to the space of vectors 
along ji("f) annihilating the 1-form dt. In jet-coordinates, the identification is expressed by the 
relation 



The push-forward of the imbedding A — ->■ ii(V n +i), restricted to the subspace A (7) C T(A) , makes 
the latter into a subbundle of ji(V {^)) . This gives rise to a fibered morphism 

Mi) h{vin)) 

(1.13a) 

viri) = v{n) 

expressed in coordinates as 

7 



dq k J- \dz A 



The kernel of the projection A (7) — ^ V(7), denoted by V(7), will be called the vertical bundle 
along 7. In fiber coordinates, V(ff) has local equation u % = 0, i.e. it coincides with the subbundle 
of A (7) formed by the totality of vectors of the form X = v A (X) (gfa")^nh) ■ 

The main interest of Proposition 11.11 stems from the study of the infinitesimal deformations. Still 
referring to pQ for details, we recall that an admissible infinitesimal deformation of an admissible 
section 7 : R — » V n +i is expressed as a section X : R — > V(7) — i.e. as a vertical vector field along 
7 — whose jet-extension ji{X) "factors" through diagram ([1. 131 a) . namely such that there exists 
a section X : R — > A (7) satisfying ji(X) = i*X. 

Under the stated circumstance the section X, viewed as a vector field along 7, is called an 
admissible infinitesimal deformation of 7. Conversely, a necessary and sufficient condition for a 
section X : R — > -A (7) to represent an admissible infinitesimal deformation of 7 is the validity of 
the relation 

(1.14a) i*X = ii(7T*(!)) 
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In coordinates, setting X = + T A (t) (gfj)- and recalling the representation (|1.13b ) . 

eq. (|1.14l a) takes the familiar form 

commonly referred to as the variational equation. 

The previous arguments point out the perfectly symmetric roles played by diagram (|1.5a ) in 
the study of the admissible evolutions and by diagram (|1.13a ) in the study of the infinitesimal 
deformations, thus enforcing the idea that the latter context is essentially a "linearized counterpart" 
of the former one. 

(ii) Let V*(j) — > R be the dual of the vertical bundle V(7) — > R, naturally isomorphic to the 
pull-back on 7 of the phase space V*(V n +i) described by the commutative diagram 

v*d) ► v*(v n+1 ) 



(1.15) 



R — ^ V n+ l 



We preserve the notation ( , ) for the pairing between V(7) and V*(j) . The elements of V*{^) 
will be called the virtual 1 -forms along 7. More generally, each element belonging to a fibered 
tensor product V(7) ®r V*{^) ®R • • • of copies of V(j) and V*(j) in any order will be called a 
virtual tensor along 7. 

Notice that, according to the stated definition, a virtual 1-form A at a point 7(t) is not a 1-form 
in the ordinary sense, but an equivalence class of 1-forms under the relation 

(1.16) A ~ A' A - A' oc {dt) j{t) 

The equivalence class determined by the 1-form dq l will be denoted by & l . With this notation, 
every section W : R — > V(7) ®r V*{^) ®r • • • is locally expressed as 

(1.17) W =w i j ...(t)(4- i 



dq l 

(iii) Given an admissible evolution 7 : R — > Vn+i , an infinitesimal control along 7 is a linear 
section h : ^(7) — > A (7) . The image ^(7) := h{V {"<()) is called the horizontal distribution along 7 
induced by h. Every section X : R — > A(fy) satisfying X(t) £ T~i( A f) Vi G R is called a horizontal 
section of -A (7) ■ 

In fiber coordinates an infinitesimal control is locally represented as 
(1.18) v A = h A (t)u l 

The associated horizontal distribution is locally spanned by the vector fields 

d_ 

dq % 



(1.19) di := h 



+h t * (± 

dq 1 )^ 1 \dz A /ry 



Every section X = X l (t) (^r) of ^(7) may be lifted to a horizontal section h(X) = X 1 di of 
-A (7). More generally, every section X = X l {t) (^?)- + T A (t) (gf?)- of A (7) may be uniquely 
decomposed into the sum of a horizontal and a vertical part, respectively described by the equations 
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(1.20a) 
(1.20b) 



V H (X) := h(n*(X)) = X*di 
V V (X) = X - V H (X) = {F A 



X*hi A ) 



d 

dz A 



(iv) A section X : R — )■ V{p/) is said to be /i-transported along 7 if the horizontal lift h(X) is an 
admissible infinitesimal deformation of 7, i.e. if it satisfies the condition i* ■ h(X) = ji(X). 
In coordinates, this amounts to the requirement 



(1.21) 



dX l 



dip 1 
dq k 



A (df_\ 
\dz A J A 



+ h k 
7 V^V 7 



X k = X k d k ip i 



By Cauchy theorem, the /i-transported sections form an n-dimensional vector space Vh, iso- 
morphic to each fiber V(j)\t. This provides a trivialization of the vector bundle V(j) — > K,, 
summarized into the identification V(j) ~ II x Vh- 

The dual space Vh* gives rise to an analogous trivialization V*(j) ~ R x Vh* ■ 

(v) The notion of /i-transport induces a R-linear absolute time derivative for vertical vector fields 
along 7. The operation may be extended to a derivation of the algebra of virtual tensor fields along 
7, commuting with contractions. 

In coordinates, introducing the temporal connection coefficients 



(1.22) 



Tk 



hi 



dz A 



we have the representation 



(1.23a) 

with 
(1.23b) 



D 
Dt 



Dt 



d_ 

dq l 



UJ J 



DW l j... ( d 



Dt 



dq i 



dW l 



dt 



+ T k i W k j ... - r/n„. + 



The algorithm is greatly simplified referring the fibers of both bundles V(j), V*{^) to 



/i-transported dual bases e( tt ) = e (l)(^r) 7 1 e^ = e 



(a) 



Setting W = W a b... e( j <8>e ( ^ ® • • • (with W a b... = W l j ... e L ■ ■ •) we have the representation 



(1.24) 



DW dW a b ... 



Dt 



dt 



Ha) 



(vi) In view of eqs. (|1.20l a.b). every infinitesimal deformation X of the section 7 admits a unique 
representation of the form X = h(X) + Y, where: 

• X = ir*X := X*(t^-)^ is a vertical field along 7, namely the unique infinitesimal defor- 
mation of 7 lifting to X ; 



• Y = VyiX) :=Y A > 9 



dz A 1 7 



is a vertical vector field along 7. 
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On account of eqs. (|1.23l a.b). the variational equation (|1.14l b) takes then the form 
written more conveniently as 



or also, in self-transported bases 

< LMb > lr=<\^l Y 

In the absence of corners, every infinitesimal deformation is therefore determined, up to initial data, 
by the knowledge of a vertical vector field along 7, through the equation 

(1.26) X a (t) = X a (t ) + j\^(^jY A dT 

1.3. Extremals. In this Section we review the stationarity conditions for a given variational prob- 
lem. Once again, the reader is referred to [1] for the full argumentation. 

(i) Let J£ € ^{A) denote a differentiable function on the manifold A, henceforth called the 
Lagrangian. Constrained calculus of variations deals with the study of the extremals of the func- 
tional I[y] := f-J^dt among all admissible sections 7 : [io>^i] — ^ V n +i connecting two given 
configurations. We shall refer to this as the control problem. 

It may be showed that, under suitable regularity assumptions, this is equivalent to a free varia- 
tional problem on the contact bundle C(A). More specifically, denoting by v : C{A) — > V n +i the 
natural projection, the control problem may be converted into to the study of the extremals of a 
suitable action functional defined on differentiable sections 7 : [io^i] C(A) whose projections 
v -7 : [tcb^i] — > Vn+i connect the given configurations. 

As a matter of fact, the contact geometry and the presence of the Lagrangian function provide 
two "natural" variational problems on C(A) . The former, of a purely geometric nature, is entirely 
determined by the Liouville 1-form (|l,10p through the action functional 

(1.27) *ra-/ s -j(M£ 

The corresponding extremals 7 : q l = q l (t), z A = z (t), pi = Pi(t) are locally described by the 
equations 

(1.28a) ^ = ^%q\z A ) 

(1.28b) f + % n - 

dip* 

(1.28c) P*^A=° 

Their role is clarified by the following observations: 

• a section 7 : [io)*i] Vn+i is admissible if and only if the functional (|1.27p admits at least 
one extremal 7 projecting onto 7 i.e. satisfying v ■ 7 = 7; 
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• the totality of extremals 7 projecting onto an admissible section 7 form a finite dimensional 
vector space £5(7), identical to the space of solutions pi = Pi(t) of eqs. ([1.28b . c) for fixed 
q l (t),z A (t). 

A section 7 is called normal if dim p(7) = 0, abnormal dim £5(7) 7^ 0; it is called locally normal if 
its restriction to any closed subinterval [to>**] — [*o>£i] is a normal arc, i.e. if and only if along any 
such subinterval eqs. Q1.28b , c) admit the one trivial solution Pi(t) = 0. 

An illustrative example of a normal evolution 7 which happen to be abnormal when restricted 
to a subinterval [tg^i] C [io>*i] may be given by means of a bump function as follows. 

Example 1.1. Let V n +l = R x £3, referred to coordinates i, g 1 ,q 2 ,q 3 , be the configuration 
space-time of an abstract system 23 subject to non-holonomic constraints. Suppose the imbedding 
A — % ii(V n +i) is locally expressed as 

q 1 = z 1 , q 2 = z 2 , q 3 = g(t) z 2 
being g(t) a C°°-function defined as g(t) := — ^sz^p e' 2 - 1 for any \t\ < 1 and g(t) := otherwise. 
Consider now the evolution 7 consisting of the arc: 

7 : g 1 = v t 2 , q 2 = v t , q 3 = vf(t) t < i ^ t\ , t < -1 , t\ > 1 

e*^ |i| < 1 

|t| ^ 1 

may be easily checked that, when restricted to the subinterval [to,— 1], equations il.28\ b.c) 
admit the family of solutions of the form 

Pi(t) = 0, p 2 (*) = 0, p 3 (t)=a, VaGR 

i/ras entailing the abnormality of the restriction 0/7 to [to,— 1]- Notwithstanding 7 is normal, 
since no solution may be found along the whole of it, but the trivial one. 

The latter variational problem is obtained lifting the Lagrangian if to a 1-form over C(A) 
according to the prescription:^ 

(1.29) : = if dt + = (Jgf - ^ dt + pi dg* := - dt + Pi d^ 

The function 34?(t, q i , z A , pi) := Piij) 1 — i? is commonly referred to as the Pontryagin Hamiltonian. 
The extremals of the functional I [7] := J- $jf are locally described by the Pontryagin equations 

dq* l A d Pl dip k d% dtf dj? 

(L30) dt = ^' Z } ' dt + ^ fc = ^' Pi d^ = d^ 

The formal analogy between eqs. (jl.3U|) and the ones arising from the solution of the original 
control problem allows to set up a relationship between the extremals of Zpy] in C(A) and those of 
the functional I[y] := L ££dt in V n+ \. 

A detailed analysis of this point can be found in [1J : for the present purposes, it is enough recalling 



/(*) == 



^To simplify things, we are not distinguishing between covariant objects in A and their pull-back in C(A) , i.e. we 
are writing tp l for (*(?/> l ), Co 1 for C{& % ) etc - 
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that every extremal 7 : [to,ii] — >• C(A) of 7(7] projects onto an extremal 7 : [to^i] ~~ * Vn+i of the 
control problem. 



(ii) A point £ € C(^4) is called regular if and only if the relation 

can be solved for z , . . . , z r in a neighborhood of (, giving rise to local expressions of the form 

(1.31) z A = z A (t,q 1 ,...,q n ,p l ,...,p n ) 

A sufficient condition for this to happen is the validity of the condition 

(1-32) ** (fl?fl?O c * 

In a neighborhood of each regular point, substituting eq. (|1.3ip into the first two eqs. (|1.30[> 
and introducing the notation 'H(t,q l ,Pi) '■= Jif(t,q l ,pi,z A (t,q l ,pi)), Pontryagin equations can be 
written in Hamiltonian form as 

1 ' J dt d Pi ' dt dq* 

In this way, the actual constrained variational problem is locally reduced to a free Hamiltonian one. 



2. The sufficient conditions 

2.1. Lagrangians "adapted" to a given extremal curve. A function / on a differentiable man- 
ifold M is said to be critical at a point x £ M if and only if its differential vanishes at x. Under the 
stated circumstance, the hessian of / at x determines a bilinear map (d 2 f) x : T X (M) x T X (M) — > H 
and has therefore a tensorial character. Similar conclusions hold if / is critical at each point of a 
submanifold N C M, in which case we write (df)w = and denote by (c? 2 /)tv the hessian of / 
along N. 

Denoting by / := |£ + ip k 6 J^{A) the symbolic time derivative of any function / = f(t, q l ) G 
j£~(V n +i), we have then the following properties: 

• if / is critical on an admissible section 7 : 1R — > V n +i , the function / is critical on the lift 
7 of 7, and satisfies = 0; 



CONSTRAINED VARIATIONAL CALCULUS: THE SECOND VARIATION (PART I) 11 

under the same assumption, for any admissible deformation X : H — > V(7) the quadratic 
form <(d 2 /) 7 ,X0X) = (J^fk) satisfies the relation^] 



7 



(2.1) ^((d 2 /) 7 ,^®^) = ({d 2 f)^X®x) 

The previous arguments are especially relevant in a variational context. For any / G J^"(V ra +i), 
the action integrals f^Jfdt and — f)dt have actually the same extremals with respect to 

fixed end-points deformations; in particular, every extremal 7 yielding a minimum for the former 
integral does the same for the latter and conversely. 

To explore the consequences of this fact, we consider once again the composite fibration 
v = -n-C, : C(A) — > V n +i described in §L3]and recall that for each normal extremal 7 : [to, h] — > Vn+X 
of the action integral f~J£dt there exists a unique extremal 7 : [ioj^i] — ^ C(A) of the functional 
f~ projecting onto 7, i.e. satisfying £ ■ 7 = 7, whence also v ■ 7 = tt ■ 7 = 7. 

In coordinates, setting 7 : (f = q % {t), z A = z A (t), pi = Pi(t), the determination of 7 relies on 
Pontryagin equations f j 1 . 30 j) . From these, it is easily seen that the evaluation of the functions Pi(t) 
is not invariant under gauge transformations =§f — > ££ — /, but depends explicitly on the choice of 
the Lagrangian. This is easily understood observing that, for Jzf' = ££ — f, eq. (|1.29p provides the 
identification 

:= Pi J + (if' + f)dt = ( Pi - + tf'dt + df 

Because of that, the extremals of the functional J- "d^i differ from those of f~ -d^f by a translation 

(2-2) Pi (t) ^Pi(t) = pi(t) - (*,«*(*)) 

along the fibres of C(A) A, as it was to be expected on account of the gauge invariance of the 
projections 7 = £ • 7 and 7 = v ■ 7. 

After these preliminaries, the idea is now to replace the original Lagrangian by a gauge equivalent 
one, satisfying the property of being critical along the section 7. 

To this end we state 



As a check we observe that the condition (df ) 7 = implies the identities 



d 2 f 

dq i dq J 



d 2 f 

dq l dz' 



d 3 f + o 3 f ^ k + d 2 f djt + d 2 f djt 

dq l dq : >dt dq l dqWq k dq l dq k dqi dqWq k dq l 



df d 2 f 



dt ydqtdqJ 

8 2 f 
dq l dq k 



+ 

dj> k 
dz A 



( d 2 f 



dip 1 

\dq i dq k J \dqi 



+ 



( d 2 f 

\dqidq k 



dq* 



d 2 f 

dz A dz l 



= 



Eq. (|2.ip is then easily obtained by direct computation, expressing the derivatives in terms of the components 
X 1 , T A through the variational equation (|1.14b ). 
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Definition 2.1. Given a normal extremal 7, a function S € ^"(V n +i) is said to be adapted to 7 
if and only if its pull-back on C(A) satisfies the condition^ 

(2-3) {dS). = (tf*)- 

By a little abuse of language, whenever a function S £ ,^{V n +i) is adapted to 7, the same 
terminology will be used to denote the Lagrangian function Jzf' = Jz? — 5. 



= 

7 



Theorem 2.1. i?g. li?. 3j) entails the relations 

(2.4) = [4^-5) 
Proof. In coordinates, a possible local solution of eq. (|2.3p is easily recognized to be 

(2.5) S (t,j) = PiWtf - j(t)) + f Stftdt 

Jo 

By linearity, the most general solution has therefore the local expression 

(2.6) S = S + f(t,q i ) 
being f{t,q l ) any function satisfying (df) y = 0. 

On account of our previous remarks, in order to complete the proof we have only to check the 
validity of eqs. (|2.4p for the function (|2.5p . To this end we evaluate 

(2-7) S = ^ (g* - «f (i)) + ft (V " ^) + 

Restricted to the curve 7, eq. (|2.7p yields back the relation (J? — So ) 7 = 0. Moreover, eqs. (|1.30|) . 
(|2.7|) provide the relations 

dq k J, dt +Pl WJ, WA 



dSn 



Pi 



d z A) \dz A ) 



which, together with ^ — So) 7 ] =0, imply the vanishing of the derivative Z(Jz? — Sq) 
VZ € T^,( t )(.A), and therefore the vanishing of [d(L — £)] » . □ 

For any S* adapted to 7, let us now replace the original Lagrangian jSf by the difference 
Jz?' = ££ — S. As already pointed out, both functions J2? and Jzf' are totally equivalent for 
the study of the extremals in V n +i , as well as for their classification. 

At the same time, by eq. (|2.2p . the extremals of the functionals f. and f~ "dy differ by a 
translation ft — > ft = ft — Jj£ along the fibers of C(.A). From this, taking eqs. (j2.5fl . ()2.6p into 

account it is easily seen that the unique extremal of f~ 1? y/ projecting onto 7 satisfies the local 
equation ft(i) = 0. 

In view of eq. (|2.6|) . the Lagrangian ££' is not unique, but is defined up to a restricted gauge 
transformation Jzf' — > jSf' + /, with (d/) 7 = 0. In any case, for each choice of 5, Theorem 12.11 
ensures the tensorial character of the hessian (t£ 2 Jz?' )& . 



5 As usual, we are not distinguishing between functions on V n +i and their pull-back on C(A) . 



CONSTRAINED VARIATIONAL CALCULUS: THE SECOND VARIATION (PART I) 



13 



The local representation of {d 2 ^£')^ is most easily expressed in non-holonomic bases. Setting 
& = (dq l - ^dt). = wj. , v A = (dz 



' — ^) 7 ' a straightforward calculation yields the result 



(2 



[d 2 ^ 



dq r dq k 



uf ® uT + 2 



dz A dq k 



Concerning eq. (j2.8[l . we observe the following facts: 



D A Qu k + 



J i 



d 2 £" 
dz A dz B 



v A ®v B 



J 1 



denoting the symmetrized tensor product. 



the components Gab '■- 



dz A dz E 



are invariant under arbitrary restricted gauge transfor- 



mations, and may therefore be evaluated arbitrarily choosing the function S within the class of 
solutions of eqs. (|2,3p . Making use of the ansatz (|2.5p we obtain the representation 



d 2 ^ 



dz A dz B 



Gab ~[ dz A dz B \~-[dz A dz B \~ Pi(t) 

written more conveniently as 
(2.9) Gab 

with K := pi(t) ^>*(i, q i , z A ) — j£f(£, q l , z A ) , henceforth called the restricted Pontryagin Hamiltonian. 
In view of the identification 



d 2 K 

dz A dz B 



d 2 K 




d 2 H 




_dz A dz B . 


lit) 


^dz A dz B . 


7(*) 



, the matrix (|2.9p is automatically non- 

• J 7(tj Lm w- J 7(tj 

singular along any regular extremal. 

• whenever detG^s 7^ 0, the hessian (|2.8p determines an infinitesimal control along 7, namely a 
linear section h : V(7) — > A{fy) , uniquely defined by the condition 

id 2 ^')., h(X)®Y] 



(2.10a) 







VI G V(j) , Y G V(j) 



In view of eqs. (|l,19p . (j2.8j) . the requirement ()2.10l a) is locally expressed by the relations 

d \ \ ( d 2 ^' 



(2.10b) 



dz A 



dq l dz A 



+ Gab hi 







7/ \ 1 / 7 

Under the assumption det Gab 7^ 0, these may be solved for the components hi B , providing the 
representation 

d 2 ^' 



hi 



whence also 
(2.11) 



di 



d_ 

dq> 



G 



d_ 

dq l 



AB 



dq l dz B 



G 



AB 



dq l dz B 



d 

dz A 



with G AB G BC = $ A c- 

On account of eq. (l2TTUb ). for any pair of vector fields U = Wdi + U A -^ , W = W*9i + W A -£j 
in A (7) , we have the relation 



(2.12) 



[d 2 ^') .,U®W) = N ij U i W j + G AB U A W B 
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with 



(2,3) Nkr := fl t .«,) - - (|^) 



7 v ^ / 7 \ ■> / 7 

_D 



The absolute time derivative along 7 induced by h will be denoted by . The expression (|1.22p 
for the temporal connection coefficients takes now the form 

All properties stated in eqs. (jl.23jl . (jl.24p will be freely used in the following. 

Unlike the components Gab , the full hessian (|2.8p and therefore also the associated infinitesimal 
control (if at all) and the corresponding time derivative are not gauge invariant, but depend on the 
specific choice of the Lagrangian ££' . 

To set things straight, when working in a local chart (U,h) containing 7 we shall stick to the 
ansatz (|2.5|) . thereby adopting the expression 

(2-15) = J2f - ^ (<f - <f (*)) - Pi {t) (V - d -£j - S>\, 

For later use, before moving on to the analysis of the second variation of the action functional, we 
eventually observe that, whenever 7 : [toj^l] ~~ * Vn+l is a normal extremal, the use of an adapted 
Lagrangian provides a canonical splitting of the tangent space T(C(A)) ■ 

As already pointed out, the definition of the (unique) extremal curve 7 : [iQjii] — > C(A) of the 
functional J- •Q^i projecting onto 7 is independent of the specific choice of ££' within the stated 
class of Lagrangians. In particular, the local representation of 7 satisfies the condition Pi(t) = 0. 

Moreover, in view of the nature of C(A) A as a vector bundle over A, each local section 
O : A — > C(A) given by pi(t, q\ z A ) = has an invariant geometrical meaning. For each a € O(A) , 
every vector Z G T a (C(A)) may therefore be split into a "horizontal" part Z h := 0*(*(Z), tangent 
to the submanifold O(A), and a "vertical" part Z v , tangent to the fiber C _1 (C( cr ))- I n coordinates, 
we have the explicit representation 

'Z,(dpi) a )( £-) , z h = z-z v 



dpi 

The previous algorithm interacts with another intrinsic attribute of the manifold C{A) , repre- 
sented by the Liouville 1-form ()1.10p . In this way, to each vector Z £ T a (C(A)) we may associate 
a 1-form in T*(C(A)) according to the prescription 

(2.16) z z v -». z v j (de) k = (z, (dpi) ff ) Q^-J de^) = (z, (d Pi 

At last, we observe that every element of T*(C(A)) generated by the process (|2.16p may be 
uniquely expressed as the pull-back of a 1-form A G T*, ^ (V n _|_i) . 

Collecting all results and recalling the definition of virtual 1-form introduced in § ll.2l we conclude 
that the lift algorithm 7—5-7 described above determines a bijective correspondence between vector 
fields Z along 7 and pairs (Z, A), in which Z = C*(Z) is a vector field along 7 = ( ■ 7, and 
A = (Z, (dpi)^) uiiy is a virtual 1-form along 7. 
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2.2. The second variation of the action functional. Let 7 : [foi^i] — > Vn+i be a normal (not 
necessarily regular) extremal of the action functional J. ££ dt. For any admissible deformation 7^ 
of 7 we denote by 7^ : q % = z A = ( A (t,^) the lift of 7^ to the velocity space A, by 

X = X*^). + T A (-£ T ). the infinitesimal deformat ion tangent to 7^ and by := J*» JSf the 
value of the action integral along each path 7^ . 



By the very definition of extremal we have then | = for all admissible deformations 7^ 

with fixed end-points. We now focus on the evaluation of the second derivative |^_ , commonly 

referred to as the second variation of the action functional at 7. 

The algorithm is greatly simplified replacing Jzf by the difference Jz?' = Jzf — 5 , being 
5 € ^(V n +i) any function adapted to 7 in the sense of Definition 12.11 By elementary calcu- 
lus we have in fact the identity 

1(0 := I (££' + S)dt = f 1 ^'{t,^(t,0,( A {t,0)dt + const. 
J% J to 

In view of this as well as of the relation (eL£?') 7 = 0, the condition ^| |^_ = is identically 
satisfied. As far as the second variation is concerned, a simple calculation yields the result 



(2.17) 



d 2 l 



d? 



£=0 J to 



X i Xi + 2 ( -^#-r ) X 1 T A + 



dq l dqi J- \dq l dz A / 

+ rArB ] dt = £ V 1 ® *} dt 



As already pointed out, the hessian (d 2 ££' ). is not a gauge invariant object. More specifically, 
according to eq. (|2.ip . under a restricted gauge transformation ££' — > J£" + /, the quadratic form 
(|2,2ip undergoes the transformation law 

(2.18) < [d 2 (J?> + /)] . , X ® X) = {(d 2 ^') v X® X) + j t {{d 2 f), n X® X) 

As it was to be expected, eq. (|2.18p shows that, unlike the integrand, the integral involved in the 
representation (|2.17p of the second variation has a gauge-invariant behaviour, at least within the 
class of fixed end-points deformations. 

As far as the second term at the right-hand-side of eq. (|2.18|) is concerned we notice that, due to 
the requirement (df) 1 = 0, the hessian (cf 2 /) 7 determines a virtual tensor field along 7, expressed 
in components as|j 

( 2 - 19 ) ( d2 f) ,= (jSL) a*® 



Conversely, for any symmetric virtual tensor C = Cij(t) Gj % <g> Cj 3 along 7, there exists at least a 
function / E ^(V n +i) defined in a neighborhood of 7 and satisfying (d/) 7 = 0, (d 2 /) 7 = C . 

The effect of the restricted gauge group on the representation (|2.17p is therefore reflected into 
the fact that the integrand at the right-hand-side is defined up to an arbitrary transformation of 



®See § l2.2l for notation and terminology. 
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the form 

(2.20) {(d 2 J?').,X®X) — ► ((d 2 J?') r X®X) + j t (C,X®X) 

being C any symmetric virtual tensor field along 7. 

The above arguments are further enhanced whenever 7 is a regular extremal. First of all, 
introducing the horizontal basis (|2.1ip associated with the hessian (d 2 j2?')^ and expressing X in 

components as X = X % d{ + Y A (-^x)^ > ec L- (12.12|) provides the identification 

(2.21) ((d 2 £").,X®X) = N kr X k X r + G AB Y A Y B 

In view of that, denoting by the absolute time derivative associated with the infinitesimal 
control (|2.1ip and recalling the representation (jl.25l a) of the variational equation in the basis 
di , (gfx)-> we conclude that, for any symmetric virtual tensor field C = C(t) along 7, there exists 
a restricted gauge transformation yielding the representation 

(2.22) {[d 2 (J?' + f)] v X®X) = ({d 2 J?') v X®X) + j t (C,X®X) = 



N - + ^m) ™ + 2 Cij (IS \ x3 yA + Gab yAyB 



with the components D * 3 expressed in terms of the ordinary derivatives and of the temporal 
connection coefficients T k = —di(ip k ) by eq. (|1.23p . 

On this basis we can state 

Theorem 2.2. Let 7 : [to^i] ~~ * V n +i be a normal extremal. Then, if the matrix Gab(^) is non 
singular at a point t* € (to, t±), there exist e > and a restricted gauge transformation ££' — > S£ ,J rf 
such that the hessian \d 2 (J£" + /)] -/^ has algebraic rank equal to r for t £ (t* — e, t* + e) . 

Proof. By continuity, there exists an interval [c,d] B t* where det Gab 7^ 0. We focus on that 
interval, and apply eq. (|2.22p to the arc j([c, d]) . Setting 



Y A :=Y A + G AB C lT (^] X> 



7 

and taking the symmetry of Cjj into account, eq. ()2.22p may be written as 
([d 2 (J?' + f)].,X®X) = 



" 1 Dt ~ v^;-v^ B A s " 



A <>B 



X l X J + G AB Y Y 



The thesis is therefore established as soon as we prove that the Riccati-type differential equation 

admits at least one symmetric solution C« = Cjj (t) in a neighborhood of t = t* . 
To this end, we set 



(2.24) M rs := G 



AB I d V \ ( d ^ 



dz A /- \dz B / 7 
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and denote by Cjp and respectively the symmetric and antisymmetric part of Cij . Due to 

the symmetry of M y and Nij , eq. (|2.23p then splits into the system 



DC, 



(S) 



(2.25 a) ^ M r 

Dt 



r (S) r (S) , r (A) r (A) 
ir sj ' ir sj 



+ Ni.. 



DC {A) 

(2.25b) — ^- + M rs Cl s) C (A) + M rs C [A) Cjf = 

Being the second equation linear and homogeneous in C- A \ by Cauchy theorem we conclude that, 
if we choose Cy symmetric at t = t* , there exists e > such that the solution of the Cauchy 
problem for eq. (|2.23p exists and is symmetric for \t — t* | < e. □ 

In view of Theorem 12.21 whenever det GAB(t*) 7^ 0, by a proper choice of the gauge around the 
point j(t*), the quadratic polynomial (|2.18p can be reduced to the canonical form 

(2.26) ( [d\X> + /)] , , X «s 1) = G AB Y*Y* = - (^g^ Y A Y B 



in a neighborhood of t*, K denoting the restricted Pontryagin Hamiltonian. 

Unfortunately, the purely local validity of eq. (|2.26p makes it unsuited to the study of the second 
variation (|2.17j) . which involves an integration over the whole interval [toil]- We shall return on 
this point in § 12.31 At present, we concentrate on the role of Theorem 12.21 in the identification of 
necessary conditions for a regular extremal 7 to yield a (relative) minimum for the action functional. 



In this connection, a preliminary result is provided by the following 

Corollary 2.1. Under the same assumptions as in Theorem \2.SX given any vertical vector field 

* dz A ) 7 



V = V A (jrx) - along 7 with compact support [a,b] C (t* — e, t* + e) , there exist a differentiate 



function g = g(t) not identically zero on [a,b] and an infinitesimal deformation X = X l (^-)^ + 
Y A ( y ^x)x j with support contained in [a,b] satisfying the relation 

Y A + G AB C r .(j^ X* = gV A 

Proof. Using the variational equation in the form (jl.25p . the required conditions are summarized 
into the pair of relations 



(2.27a) 



BX l 
Dt 



\dz B J, 

(2.27b) X\a) = X\b) = 

For any choice of g{t) , eq. (|2.27|) is a first order linear differential equation for the unknowns 
X l (t); integrating it with initial data X % (a) = yields the solution 

x\t) = w\(t) [(w-^K (?j^gV A dt 

being W l k the Wronskian of the equation. In order to ensure X l (b) = it is then sufficient to choose 
g(t) within the (infinite-dimensional) vector space of differentiable functions over (t* — e, t* + e) 
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satisfying the conditions 

j\w-^ r (^y 9 V A dt = 0, k = l...n 

□ 



Corollary 2.2. A necessary condition for a normal extremal 7 : [io)*i] ~~ * Ki+1 to yield a minimum 
for the action functional is that the matrix GAB{t) be positive semidefinite at all t £ [£o,ti]- 

Proof. Suppose that Gab is not positive semidefinite at some t* £ [ioj*i]- Depending on the value 
of det GAB(t*) we have then two possible alternatives: 

i) if det Gab (t*) 7^ , on account of Theorem 12.21 there exist a restricted gauge transformation 
Jzf' — > J£' + / such that a representation like ()2.26|) holds in a neighborhood (t* — e, t* + e) . 

Then, given any vertical vector field V with support contained in (t* — e, t* + e) and satisfying 
Gab y A y B < (for instance, the eigenvector corresponding to the negative eigenvalue of Gab in 
(t* — e, t* + e), multiplied by a suitable function with compact support), Corollary 12.11 implies the 
existence of at least one infinitesimal deformation X satisfying 



d 2 ! 
d£ 2 



5=0 J to 



f 1 ([d\j?' + f)] X®X)dt = f g 2 G AB V A V B dt < 

J to J a 



Therefore, 7 does not provide a minimum for the action functional, 
ii) if det Gab (O = 0, choose e > in such a way that 

• — e is not a root of the secular equation det (Gab — A#ab) = 0; 

• at least one root of the secular equation is smaller than —e. 

Let j$ G <^{A) be a differentiable function globally defined on A and having local expression 
^# = £ Sab (z — z A (t))(z B — z B (t)) in a neighborhood U of the point j(t*) , where, as usually, we 
are writing z A (t) for z A (j(t)). Also, let [c,d] 9 t* be a closed interval, satisfying j([c, d]) C £7. 
Setting Jzf* := «Sf' + one can then easily verify the properties: 

a) the section 7 : [c, d] — > V n +i is a normal extremal for the action integral J~ L*dt; 

b) the matrix ( )~(t*) = Gab+z &AB is both non singular and non positive (semi)-definite. 

In view of a) and b), the analysis developed in point i) ensures the existence of at least one 
infinitesimal deformation X having support in [a, b] C [c,d] and satisfying 

((d 2 ^*)^,X ®X)dt < 
On the other hand, by construction, this implies also 

J ^(d 2 i?')_ , X® X^jdt = J ((d 2 ^*), , X®X^dt + 

- £ J 5 A B^dz A , X^ (dz B , X^dt < J <J(d 2 if*). , X®X^dt < 
once again proving that 7 does not yield a minimum for the action functional. □ 
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2.3. Regular extremals: the matrix Riccati equation. From now onward we shall concentrate 
on the class of regular normal extremals. The role of regularity in the solution of the Pontryagin 
equations (|1.3U|) — more specifically, in the conversion of these into a system of ordinary differential 
equations in Hamiltonian form for the unknowns q l (t),pi(t) — is well-known in the literature^ 
However, in the present context, regularity is meant as a merely attribute of a given extremal 
7 : [io>£i] ~~ * Ht+ii ensuring the existence of an expression of the form (|2.26p in a neighborhood of 
each t* £ [to, *i] • 

On the other hand, as already pointed out, the purely local validity of eq. (|2.26p is of little help 
in the evaluation of the second variation (|2.17p . In this connection, it ought to be established to 
what extent eq. (|2.26p may be converted into a global result, valid over the whole interval [io,ii]- 
On account of eqs. ()2.23|) . (|2.24p . this entails to analyze the interval of existence of the solutions of 
the Riccati-like differential equation^ 

DC- 

(2.28) - M rs C ir C sj + N tj = 

The main difficulty comes, of course, from its non-linearity. To overcome this aspect, we intro- 
duce two auxiliary virtual tensors Eij(t) and K' l j(t) along 7, satisfying the transport laws 



DK 
~~Dt 

DE. 



(2.29a) — ^ = M ir E rj 

Dt 3 

Dl " j 



(2.29b) ^L = Nir K r 

On any interval (a, b) on which det K l j(t) 7^ 0, the (generally non symmetric) tensor 

(2.30) Qj =-E ir (K- 1 ) r j 

is thus well-defined and fulfils the relation 

DE^ _ DQr DK% 
Dt ~ Dt p+ Uir Dt 

Substituting from eqs. (1239k L ([Z3DD and multiplying by (K^) p 'j the latter may be written in the 
form 

Tin. . nn. . 

mS!) - N t] = + C ir M rs E sp (K- 1 )^ = -^L - C ir M rs C sj 

formally identical to eq. (|2.28p . 

Needless to say, the symmetry property Cij = Cji is also needed in order for the tensor (|2,30p 
to represent the hessian of a function / along 7. In view of eq. ()2.28f ). the antisymmetric part of 
Cij obeys a linear homogeneous system of the form (|2.25b ) and therefore an argument similar to 
the one exploited in the proof of Theorem 12.21 shows that this aspect relies entirely on the choice 
of the initial data. Once again, by Cauchy theorem we conclude that if Cij is symmetric at t = to 
(as it happens e.g. choosing Ei r (to) = 0, K r j(to) = 5 r j), it retains its symmetric character up to 
the first value t* > t (if any) at which det K r j(t*) = 0. 



7 See e.g. [J. 

^The regularity assumption is once again crucial in ensuring the global character of the absolute time derivative 
j^r induced by the hessian (d 2 ^f')^ along 7. 
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Remark 2.1. The analysis of eqs. {2.28\) . \2. 29\ a. b) is considerably simplified 

referring the virtual tensor algebra along 7 to an h -transported basis { e , e ( a ) } and recalling 

DT a 

that, in this way, the components — D £ " of the absolute time derivative of a field T coincide with 

the ordinary derivatives dT d ^"' (§ \2.2l eq. {1.2$ ) • Eq. {2. 28\) reduces then to the ordinary matrix 
Riccati equation 

(2.31) ^ - M rs C ar C sb + N ab = 

while eqs. {2.29a . b) take the simpler form 

(2.32 a) ^ = M ac E cb 

(2.32b) ^ = N ac K\ 

Collecting all results, we can now state 

Theorem 2.3 (Sufficient conditions). Let 7 : [to,ii] — > V n +i be a normal extremal for the action 
functional. Also, let ^ = piip 1 — Jz? denote the Pontryagin Hamiltonian associated with the given 
Lagrangian. Then, if the matrix 

Gab (£) := ' S ' K ^ ' 



dz A dz B J ~ \dz A dz B J - 

is positive definite at each t £ [io>ii] an d if the system {2.29a . b) admits at least one solution 
Eij(t) , K l j(t) fulfilling the conditions 

• E{ r (if -1 ) j symmetric, 

• det K l j ^0 everywhere on [io>il]j 

the section 7 yields a weak local minimum of the action functional. 

Proof. The stated assumptions imply both the regularity of the extremal 7 and the existence of a 
global solution of the Riccati-like equation ()2.28|) along 7 , thus ensuring the validity of an expression 
like eq. f)2.28|) on the whole interval [io)*i]- On account of eq. (|2.17p . if the matrix Gab is positive 
definite on [to,ti], this provides the evaluation 



d 2 l 



d? 



[ 1 ([d 2 (_§?' + /)] . , X®x)dt = I 1 G AB Y A Y B dt > 



5=0 Jt Jt 



for every non-zero admissible deformation X : [ioj^l] ~> D 

A deeper insight into the meaning of the condition det K l j 7^ is provided by the study of the 
Jacobi vector fields, reviewed and adapted to the present geometrical context in §[3j 



3. The necessary conditions 

3.1. Jacobi fields. Given a regular normal extremal 7 : [£o,ii] — > V n +i of the action functional 
L Jzf dt, we now focus our attention on the (unique) extremal 7 : [io,£i] — > C(A) of the functional 
J* projecting onto 7. Also, as usual, we preserve the notation v = n ■ £ : C(A) — > V n +\ for the 
composite fibration C(A) — > A — > V n +i- 
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Let us now introduce a special class of deformations 7? of 7 in which every section 
7£ : [to>*l] — ► C(A) is itself an extremal of f~ $3?'. 

In this way, the 1-parameter family 7^ := v ■ 75 is a deformation of the original section 7, 
consisting of extremals of the functional f~ JC dt. At this stage, we do not impose any restriction 
on the behaviour of the end-points 7g(to)> 7g(^i)- I n coordinates, setting 

(3.1) 7? : j=<P%Z), z A = ( A (t,0, Pi=Pi(t,0 

our assumptions are summarized into the request that, for each value of £, the functions at the 
right -hand-side of eqs. f)3. 1 1) satisfy Pontryagin equations 

(3.2a) ^ = ^(t, V \C A ) 

(3 - 2b) -m + W pk = W 

(3 ' 2c) Pi d^=a^ 

As a check of inner consistency it is worth observing that, in view of the condition (dJf')z = 0, 
eqs. (|3.2b . c) and the normality of 7 yield back the relation pi(t, 0) = 0. 

Strictly associated with 7* is a corresponding infinitesimal deformation, locally expressed as 

-(ft— (fl,- (ft 

Taking eqs. (|3.2p and the relation pi(t, 0) = into account, it is easily seen that the components 
satisfy the following system of differential-algebraic equations 



(3 - 4b) iF + 7Tk \w)i~ wa? ), + \m? 

(3 - 4c) Vi 8^ = Kd&d?)* + {dz^)- T 

Given any vector field X satisfying eqs. (|3.4p . its push-forward v*X will be called a Jacobi field 
along 7. By definition, a Jacobi field X = X i [^-x) is therefore the infinitesimal deformation 
tangent to a finite deformation 7^ consisting of a 1-parameter family of extremals of the functional 

Remark 3.1. The resemblance between eqs. flff.^p and Pontryagin ones il.30\) sticks out a mile. 
This aspect can be made explicit by replacing the imbedding A1.5\ a.b) with its linearized counterpart 
\1.13\ a.b). namely regarding the vector bundle V(7) as the configuration space-time of an abstract 
system 33' and the bundle A{^) — > V{^) as the associated space of admissible velocities. In this way, 
the admissible evolutions of 23' are easily seen to be in bijective correspondence with the admissible 
infinitesimal deformations of 7 . 
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Referring V{j) to coordinates t,u l and A(j) to coordinates t,u l ,v A according to the prescrip- 
tions \l.lS\ a.b). the imbedding i* : A(fj) —tjiiVfr)) is locally expressed by eq. hi. 13b ). synthetically 
written as 



(3.5) 



u 



dtp 1 
dq k 



u k + 



dip 1 
dz A 



7 \ — / 7 

To complete the picture, we adopt the quadratic form 



(3.6) 



£(X) := ({d 2 J?') r X®X 



as a Lagrangian on A (7), and denote by 3 the functional assigning to each admissible section 
X : [to,t\] — > V(7) the action integral 3[X] := j-^Qdt. 

In this way, for any finite deformation 7^ of 7 tangent to X , eq. {2.1 r ty provides the identification 

^ 1X1 = \ ^|?| ?=0 - 

In coordinates, eq. h3. 6\) takes the explicit form 



£(t,u\v J 



dq l dqi 



u l w> + 2 



dq l dz J 



u i v A + 



3 2 S£' 
dz A dz B 



v A v B 



From the latter it may be easily verified that the equations \3.4\) involved in the definition of the 
Jacobi fields are formally identical to the Pontryagin equations for the determination of the extremals 
of the functional 3 subject to the constraints |X 



Rather than pursuing the analogy outlined in the previous Remark, we work directly on the 
system (|3.4p . Recalling the discussion at the end of ^2.1| the field X may be decomposed into the 
pair 

in which X is a vector field along 7 while A is a virtual 1-form along 7. By a little abuse of 
language, this will be called a Jacobi pair belonging to X = v*X. 

Under the (crucial) hypothesis of regularity of 7, we now make use of the infinitesimal control 
h : V(7) — > A(fy) induced by the Lagrangian Jz?' to express the field X in terms of the Jacobi 
field X and of a vertical vector Y in the form X = h(X) + Y = X l di + Y A -^^ . On account of 
eqs. (|2.1ip . this implies the relation 

(3.7) T A = Y A - G AB (4^r) 



dq l dz B 



1 



Together with the identification Gab = \ q z aq z b ) 7 > the latter allows to cast eq. (|3.4b ) into the form 

(3.8) »<^)r eju,y " -* yA = gAB9 '(^ 

From this, substituting into equations ()3.4l b. c), recalling the definitions (|2.13j) . (|2.24p of the 
tensors Nij , M !J and expressing everything in terms of the absolute time derivative, we eventually 
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obtain the following system of differential equations for the unknowns X' 1 (t) , fc j (t) : 



2.3 



(3.9a) 



(3.9b) 



DX l 
~Dt 

Dt 



G 



AB 



dip 1 \ ( dtpi 



^dz A 

d 2 ££' 
dq l dqi 



8. 



,B 



7T, 



G- 



dq l dz A J^\dqidz B 



ab( d 2 J?> 



X 3 



Nij X 3 



As the attentive reader will have noticed, these are formally identical to the linearized form 
of Riccati equation (|2.29p : a result that will prove to be fundamental in the sequel in order to 
determine the necessary and sufficient conditions for a weak local minimum. 

Remark 3.2. Keeping in line with Remark \3.1l if the virtual tensor algebra along 7 is referred 
to an h-transported basis {e^ a \e( a )} ; the set of 2n differential equations \3. 9\) are written in the 
form 



(3.10) 



dX a 
~dT 



M ab fr b 



dir a 
dt 



N ab X b 



Once again, these are easily seen to represent the Hamilton equations for the Hamiltonian func- 
tion 



(3.11) 



fj(t,X a ,ir b ) = 7f a tf° - £ 



1 M ab TT a TT b - \N ab X a X b 



2 2 

The relationship between Jacobi fields and the second variation is made evident by the following 



Proposition 3.1. Given a Jacobi pair (X,X), any infinitesimal deformation 

<■ dz A ) 7 



Z = Z % {pi)~ + Z A (-r^x)~ satisfies the relation 



(3.12) Ud 2 %')^X®Z 



7' / dt 

Proof. Eqs. (jl.25l a) and (|3.9l b) provide the evaluation 

^ = = N ijX ^ + *Mz A = N ijX ^ + GAB Y°Z* 

dt Dt oz A 

The thesis immediately follows in the light of eq. (|2.12p . □ 

Remark 3.3. Hitherto, our treatment of Jacobi fields has uniquely involved the adapted Lagrangian 
Jzf' . This choice was suggested both by consistency with the previous analysis and also by the 
simplified calculations. However, it is not at all necessary in order to cover the subject. The same 
considerations can be obviously done sticking on the original Lagrangian function. In this way, 
equations (|3.4p are directly written in terms of the Pontryagin Hamiltonian J^f , with the quantities 
7fj replaced by 7Tj := (^)^_ and related to the previous ones by the relation 

fdpi(t,0\ 8( . „ dS\ _ d 2 S . 

The argument is almost identical to the one developed so far and will be omitted. 
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3.2. Conjugate points. The arguments of § 13.11 come helpful for establishing the necessary con- 
ditions for a given extremal to yield a local minimum for the action functional. 

Definition 3.1 (Conjugate point). A point 7(7"), r G (to,ti\, along a given extremal curve 
7: [toj^l] ~~ > Vn+i is said to be conjugate to j(to) if and only if there exists a non-zero Jacobi 
field X : [t ,ti] ->• 1/(7) such that X(t ) = X(t) = 0. 

The search for conjugate points can be performed by looking for a solution of equations (|3.9[) with 
X l (to) = and vfj(to) varying amongst all the possible values in K™ . Because of the linearity 
of equations (|3.9p . their solution will depend on the initial data through a set of time-dependent 
matrices in the form 

(3.13a) X\t) = srf){t,t G )Xi{t ) + t ) in,- (to) 

(3.13b) 7Vi(t) = ^ > ij (t,t )X^t Q ) + ®i j {t,to) TTjito) 

with £/)(t Q ,t ) = ^/(t ,io) = 5) , ^(t ,t ) = %j(t ,t ) = 0. 

Hence, by Definition 13.11 a point 7(7-) is conjugate to 7(^0) if and only if the map 
V* (7)|t —> ^(7)1* locally described by 

(3.14) X\t) = ^(t,t )7rj(t ) 

is injective for t G (to,r) but not for t = r, i.e. if and only if the associated matrix satisfies 
det B^(t,t ) + for t Q < t < r, and det /3 ij '(r,t ) = 0. 

The link between conjugate points and the analysis of the second variation is clarified by the 
following generalization of a classical result of Bliss (|25J: 

Theorem 3.1. Consider an extremal curve 7 : [io^i] ~~ * V n +i and suppose there exists a value 
t G (to,ti) such that the point j(t) is conjugate to 7(^0)- Then, the quadratic form 

(3.15) f 1 ((d 2 ^'),,X^x)dt 



J to 

is indefinite within the class of infinitesimal deformations X vanishing at the end-points. 
Proof. Let us define a symmetric bilinear functional (<i 2 Z)^ over ^(7) as 

Ud 2 X). , V®w) := t /{d 2 J?')„, V®w\ dt 

J tQ 

for any V, W in ^(7). Then, by a well-known result in the theory of quadratic forms, the thesis 
is proved as soon as we show that, in the presence of a point 7(7-) conjugated to 7(^0) > the kernel^ 
of (d 2 l) - does not coincide with the locus of zeroes of its associated quadratic form. 

Under the stated hypothesis, there exists a non-zero Jacobi field J G V{^) such that 
J (to) = J(t) = 0. By means of this, we now define a continuous infinitesimal deformation vanishing 



^We recall that the kernel of the bilinear functional (d 2 T) ^ : A (7) x A (7) — 5* R is defined as the set ker(cPl) . 
{V\V e A(j),((d 2 l) . , V ® W) = VW g A(j)}. 
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at the end-points X : [io>*i] ^(7) as follows: 

J(t) t ^t^T 



X(t) :-- 



T € t < U 



Then, preserving the notation X for the lift of an infinitesimal deformation X and denoting by 
( J, A) a Jacobi pair belonging to J, in view of equation (|3.12p we have 

d 2 l).,X®x) = f Ud 2 jSf"). , J ® j) dt = 



to 



Tti J 1 



= 

to 



At the same time, if W is any infinitesimal deformation of 7 vanishing at the end-points and such 
that W(t) 7^ 0, we have also 

{d 2 l), , x®w) = jT ({d 2 ^'), , 3®w)dt = 

Since, by hypothesis, J(t) 7^ for every t E (to,r), the uniqueness of the solution of the "time- 
reversed" Cauchy problem (|3.9|) in 7(7") implies that 7fi(r) 7^ for at least one value of the index i. 
Therefore 

' U 2 z). , x®w) ^ 



7 

showing that X does not belong to the kernel of (<i 2 Z)^ . 

This completes the proof, as the indefinite character of the quadratic form (d 2 Z)^ is made 

apparent simply by looking at its action on any linear combination of the form aX + W , for any 
a G R. □ 



As an immediate consequence of Theorem 13,11 we have the following 

Proposition 3.2 (Necessary conditions). Suppose the extremal closed arc 7 : [to>£l] V n +i 
is a (local) minimum of the functional I[y] = Lj£dt. Then, for every r G (to,t\), there cannot 
be any point 7(7") conjugate to 7(to)- 



4. Necessary and sufficient conditions 

4.1. Riccati and Jacobi equations. So far we have separately proved a sufficient and a necessary 
condition for a given regular normal extremal 7 to be a minimum; we shall now merge them together 
into a necessary and sufficient one. In the event, we will prove that, whenever no conjugate point is 
present, the solutions of Jacobi equation can be used to build a global solution of Riccati equation, 
thereby allowing to apply Theorem 12.31 to the whole interval [io^i]- 

However, in order to do so, we shall need to strengthen the hypothesis of normality of 7 by 
requiring the latter to be locally normal. 
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First of all, recalling that the closed arc 7 : [to, t\] — > V n +i is the restriction to the closed interval 
[io>*i] °f an admissible section defined on some open neighborhood (a, b) D [£o>*i]> we consider 
the following technical argument 

Lemma 4.1. Let 7: [toi^i] Ki+l be a locally normal extremal and suppose the matrix Gab is 
non-singular at each t E [ioj *i] - If along 7 there is no point conjugate to 7(io)> then there exists 
a t* > t\ such that the absence of conjugate points may be extended over a wider interval \t$,t*\. 

Proof. Consider the family of Jacobi pairs (Xr^ , A(^) , k = l,...,n, obtained as solutions of 
equations (j3.9j) with initial data 

X \k)^) = ' *i(fc)(*o) = $ik 

The non-existence of conjugate points along 7 is easily seen to be equivalent to the condition 
det (X* (k) (t)) ^0 for all iE(t ,*i]- 

If that is not the case, there would be some r E (to,t\] at which the homogenous system 
a k X^Jt) = would admit a non-null solution a 1 ,... ,a n . The fields X := a k X^, A = a k 

would then constitute a Jacobi pair satisfying the conditions A (to) ^ 0, X(io) = -X'(r) = 0. On the 
other hand, X cannot be identically zero over the whole interval [to> r ] : if it were so, the 1-form 
A would satisfy the equations 



» M J 7f 1 =0 7f,- = 



Dt A, " J " 

= 

7 



V to ^ t ^ r 



L>t 



contradicting the local normality of 7. 

To sum up, X would be a non-zero Jacobi vector field vanishing at both to and r, which clashes 
with the assumption of non-existence of conjugate points along 7. 

By continuity, this implies the condition det (X\^ (t)) 7^ for all t E (toi**] with t* E (ti,&) 
sufficiently close to ti. The absence of conjugate points holds therefore in a wider interval [iojt*]- 

□ 



Coming back to the formulation of the necessary and sufficient conditions for minimality, we are 
now able to state the following 

Proposition 4.1. Let 7 : [io>*i] ~~ > V n +i be a locally normal extremal and suppose the matrix Gab 
is non-singular at each t E [to>^i]- If no pair of conjugate points exists on 7, the Riccati equation 
\2. 28\) admits a symmetric solution throughout the interval [io>*i]- 

Proof. As usual, we regard 7 as the restriction of an admissible section defined on an open interval 
(a, b) D [io ?£].]• Let t* E (ti,b) and consider a family of solutions (Xr^, Xr k \) of equations (|3.9p . 
obtained imposing the initial conditions XJ k Jt*) = and (7f)i(jfc)(i*) = &ik- 

In view of Lemma 14.11 whenever t* is chosen sufficiently close to t\, the absence of conjugate 
points implies the requirement det (X^Jt)) 7^ for all t E [to,**)- 
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A comparison between the Hamiltonian system (|3.9p and the linearization (|2.29p of Riccati 
equation shows that we can now assume the identifications 

(4.1) K^it) := X\ 3) {t) , Eij(t) := n i{j) (t) 

As a consequence, the matrix K l j(t) is non-singular everywhere on [to,t*) and therefore, as we 
have seen in § 12. 3| the tensor Cy = Ei r [K" 1 ) j represents a solution of Riccati equation (|2.28|) 
all over the interval [to,t*) D [to 5 *i] • 

In order to complete the proof, we now only need to show that the stated CL- is also a symmetric 
tensor. To this end, we first observe that the matrix R %p := K % j (E~ i y p is perfectly meaningful in 
a neighborhood (t* — 6, t*\ and satisfies the relations 

R ip {t*) = , R rp C pq = K { j (E~ 1 y p E w [K^Yq = S% Vt < t* 

The matrix R ip is therefore symmetric at t = t* . Moreover, on account of equations (|3.9p . it 
satisfies the equation 

EH! = £Ei( E -yp + jc D ( E ~ 1 Y P = W r _ R u NlkR k P 

Dt Dt y ' 3 Dt 

which is again of the Riccati-type ()2.28p . with the roles of the matrices M lJ , Nij interchanged. 
Exactly as in Theorem 12.21 this establishes the symmetry of R tp in a neighborhood of t = t* . 

For each t € (t* — 5,t*) the matrix Cij(t) = 1 is therefore symmetric. Once again, on 

account of the linearity of equation (|2.25l b). we conclude that Cy(t) is symmetric over the whole 
interval [to,t*) D [io^i]- ^ 

Collecting all the above arguments, we are now able to state the following 

Theorem 4.1 (Necessary and sufficient conditions). Suppose the closed arc 7: [io>^i] ~~ ^ Vn+i 
is a locally normal extremal of the functional I[y] = f-JZ'dt with respect to the class of deformations 
vanishing at the end-points and let 7 be its ( unique ) lift to C(A) solving Pontryagin equations 
Denote by Jif(t,q l ,z A ,pi) = pi^ l (t,q l ,z A ) — Jjf (t, q" 1 , z A ) the Pontryagin Hamiltonian associated 
with the given Lagrangian and let 



G AB {t) :-- 



dz A d 



7 

The section 7 is a minimum of the functional I [7] if and only if, for any t € [io>^iL 
the matrix Gab is positive definite and there is no point conjugate to 7(io)- 

The proof should, at this time, be quite straightforward and is left to the reader. 
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